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Abstract 

We measured two moments of the lepton momentum spectrum in B — > Xiv, where £ = e oi fi, for 
Pi > 1.5 GeV/c. From these we derive the Heavy Quark Expansion (HQE) parameters A(MS) = 
0.39±0m\stat±0-0Q\sys±0.12\th GeVand Ai = -0.25±0.02\stat±0m\sys±0.U\th GeV^, through 
order in the non-perturbative expansion and Poa^ in the perturbative expansion. The 

theoretical expression needed to extract \Vcb\ from the measured semileptonic width is evaluated 
using these HQE parameters. Combined with the world average of the semileptonic width, we find 
\Vcb\ = (40.8 ± 0.5|r^, ± ^■M{x-i_.A)^^p ^ 0-9|t/i) x 10~^. Finally, the short range 6-quark mass m^^ is 
evaluated and found to be 4.82 ± 0.07|exp ± 0.11\th GeV/cP. 



*Submitted to the 31®* International Conference on High Energy Physics, July 2002, Amsterdam 
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I. INTRODUCTION 



Experimental data on inclusive B meson semileptonic widths have reached great preci- 
sion [0. In order to use these data to determine the Cabibbo-Kobayashi-Maskawa (CKM) 
parameter \Vcb\, we need a theoretical evaluation of the hadronic effects for this process. 
We can make use of theoretical expressions for the differential rate dTgi/dE^ in terms of 
\Vcb\- The heavy quark expansion [Q, |, |, |^ is a QCD-based approach that casts perturba- 
tive and non-perturbative corrections to the partonic width as power series expansions. An 
underlying assumption of this approach is quark-hadron duality. 

It is important to assess the uncertainty introduced by possible duality violations quan- 
titatively, in order to achieve a full understanding of the theoretical errors and be able to 
ascertain the true uncertainty on \Vcb\- 

The theoretical expression for the inclusive semileptonic width for B Xiu {£ = fi or 
e) through 0(1/M|) in the non-perturbative expansion and /3o(— )^ in the perturbative one 
is given by ||, g 



r,; = ^ ' ^ 0.3689 

1927r^ 



1.54^-1.43/3of-'' 

TT V TT 




— (l- 0.87^) - 0.946 f— V - 3.185^ + 0.02^ 

A \ AiA AoA pi oo 

' 3.28^^ + 10.47^^ -6.153-^ + 7.482- 



Mb; "Mb'^ ■ Mb^ ■ Mb^ ■ Mb^ 

-7.4^^ + 1.491^^ - 10.41^^ - 7.482^^ + O I ^-^ 

Mb^ Mb' Mb' Mb' \Mb\ 



The parameter Ai Q is related to the average kinetic energy of the h quark inside the B 
meson: 



>^^ = l^^{B{v)\hSDYK\B{v)) . (2) 

The parameter A2 0, |^ is the expectation value of the leading operator that breaks the 
heavy quark spin symmetry. It is formally defined as 

^^ = W^ (B{v)\hJ- ■ a^^G,MB{v)) , (3) 



2Mb \ " " 2 

where hy is the heavy quark field, and \B{v)) is the B meson state. It is determined from 
the B* — B mass difference to be 0.128 ± 0.010 GeV^. A is related to the 6-quark pole mass 
nib 0, ^ through the expression: 

rrtfe = Mb - A + (4) 
2mb 

where Mb is the spin-averaged -B*-*^ mass {Mb = 5.313 GeV/c^). A similar relationship 
holds between the c-quark mass rric and the spin-averaged charm meson mass (Md = 1.975 
GeV/c2). 

The shape of the lepton energy spectrum in the decays B — > Xiu can be used to measure 
the HQE parameters Ai and A, through its energy moments, which are also predicted in the 
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Heavy Quark Expansion. We choose to study truncated moments of the lepton spectrum, 
with a lepton momentum cut of pe > 1.5 GeV/c. This choice decreases the sensitivity of our 
measurement to the secondary leptons from the cascade decays (6 — c — > s/diu). 
The quantities that we focus on, originally suggested by Gremm et al. [0,^ are 

_ J,^,Ei{dTsi/dEi)dEi 

' j,^,{drjdEi)dEi • ^ ^ 

The theoretical expressions for these moments R*'^ 0, H are evaluated by integrating over 
the lepton energy in the decay b — > c£u for the dominant Fc component. In addition, the 
small contribution coming from charmless semileptonic decays b —>■ uiu is included in the 
theoretical expression by adding the integral over dTu/dEi where the contribution to the 
moments is scaled by J^j HI- 

We determine these two moments from the measured lepton spectrum in B ^ Xiu 
and insert them in the theoretical expressions to determine the two parameters Ai and A. 
Previous experimental determinations of A and Ai have been obtained by studying the E^ 
spectrum in 6 — >■ 57 and the first moment of the mass Mx of the hadronic system recoiling 
against the iu pair in i? ^ Xiu decays |T0|. We will summarize the present knowledge on 



these parameters and the implications for |V^b|. 

In recent years, increasing attention has been focused on "short-range masses," preferred 
by some authors as they are not affected by renormalon ambiguities [^. In particular. 



theoretical determinations of the so-called fe-quark mass in the 15* renormalization scheme, 
mP, based on extractions from T(1S) resonance data are quoted with impressive theoretical 
accuracy []T2|. Using the formalism developed by C. Bauer and M. Trott we have used 



semileptonic energy moments to extract m 



IS 



II. EXPERIMENTAL METHOD 

The data sample used in this study was collected with the CLEO II detector [^] at the 
Cornell e'^e" collider (CESR). It consists of an integrated luminosity of 3.14 fb~^ at the T(4S) 
energy, corresponding to a sample of 3.3 x 10^ BB events. The continuum background is 
studied with a sample of 1.61 fb~^ collected at an energy about 60 MeV below the resonance. 

We measure the momentum spectrum of electrons and muons with a minimum momentum 
requirement of 1.5 GeV/c. This momentum cut serves two purposes. First, it ensures 
that we can take advantage of good particle identification criteria for both lepton species. 
Having both lepton species allows us to check crucial systematic effects by examining the 
fi/e ratio. In addition, in this range the inclusive spectra are dominated by the direct 
b ciu semileptonic decay, rather than by secondary leptons produced in the decay chain 
b —>■ c —>■ siu. 

Electrons are identified with a likelihood method that includes several discriminating 
variables, most importantly the ratio E/p of the energy deposited in the electromagnetic 
calorimeter to the reconstructed momentum, and the specific ionization in the central drift 



Our notation is different than what is used in Rcf. ||7|, where Rq is first introduced as R2. 
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chamber. Muon candidates are required to penetrate at least five nuclear interaction lengths 
of absorber material. We use the central part of the detector (| cos6'| < 0.71 for electrons 
and I cos^l < 0.61 for muons). 

The overall efficiency is the product of three factors: the reconstruction efficiency, includ- 
ing event selection criteria and acceptance corrections; the tracking efficiency; and the /z or e 
identification efficiency. The first two factors are close to 1 and estimated from Monte Carlo 
simulations and checked with data, whereas the lepton identification efficiencies are studied 
with data: radiative /i-pair events for the fi efficiency and radiative Bhabha electron tracks 
embedded in hadronic events for the e efficiency. The e identification efficiency is rather 
constant in our momentum range and equal to 93.8 ± 2.6%. The fi momentum threshold is 
near our low momentum cut, but the efficiency rises to a plateau of about 95% above 2.0 
GeV/c. 

Figure |I| shows the raw yields for electrons (top) and muons (bottom) from the T(4S) and 
the scaled continuum background. The scaling factor for the continuum sample is determined 
by the ratio of integrated luminosities and continuum cross sections and is 1.930 ± 0.013. 
This scale factor has been determined independently using high statistics samples of tracks 
with momenta higher than the kinematic limit for S-meson decay products. This control 
sample has been used to determine the systematic errors in the measurement. 

The raw yields include hadrons misidentified as leptons (fakes). This contribution is 
determined from data as follows. Fake rates are determined from tagged samples: tt's 
from Kg — >■ 7r+7r~, K^s from D*~^ D^7r~^,D^ — > K~7r^, and p's from A — >■ pir^. Using 
the measured fake rates and the observed spectrum of hadronic tracks, we compute a fake 
contribution, which is quite small for both electrons and muons. 

Several sources of real leptons must be corrected for. Leptons from J/ip decays are 
vetoed by combining a candidate with another lepton of the same type and opposite sign 
and removing it if their invariant mass is within 3a of the known J/ip mass. Corrections are 
applied for veto inefficiency and for signal leakage into the veto window. Similar procedures 
are applied to electrons and positrons coming from Dalitz decays and from 7 conversions. 

Finally, we subtract leptons coming from ip' decays or the secondary decays b ^ c ^ 
siu and B r Ivv using Monte Carlo studies. Figures ^ and |^ show the individual 
estimated background contributions to our sample. Note that all of the backgrounds are 
small compared to the signal. 

Our goal is a precise determination of the shape of the lepton momentum spectrum, so 
corrections for the distortion introduced by electroweak radiative effects are important. We 



use the prescription developed by D. Atwood and W. Marciano [|14[. This procedure incorpo- 



rates leading-log and short- distance loop corrections, and sums soft- virtual and real-photon 
corrections to all orders. It does not incorporate hard-photon bremsstrahlung, which mainly 
modifies the low energy portion of the electron spectrum, and is not used in our analysis. An 
independent method of studying QED radiative corrections in semileptonic decays, based on 
the simulation package PHOTOS [0, has been used to obtain an independent assessment 
of the corrections. The difference between the two methods is used to obtain the systematic 
error of this correction. 

Finally, we use a Monte Carlo sample of signal events to derive an unfolding matrix to 
boost the spectra into the 5-meson rest frame (B mesons typically have a momentum of 
^ 300 MeV/c in the laboratory frame). Figure ^ shows the resulting electron and muon 
spectra. 

Our first goal is to extract the truncated moments Rq and Ri defined in Eqs. |^ and |^. 
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Using the measured spectra, we evaluate the relevant integrals and obtain the results shown 
in Table |I|, where the first error is statistical and the second is systematic in each quoted 
number. 





Kg 






0.6184 ±0.0016 ±0.0017 


1.7817 ± 0.0008 ± 0.0010 


/^^ 


0.6189 ± 0.0023 ± 0.0020 


1.7802 ±0.0011 ±0.0011 


Combined 


0.6187 ±0.0014 ±0.0016 


1.7810 ±0.0007 ±0.0009 



TABLE I: Measured truncated lepton moments for e and n , and for the sum. 

The dominant uncertainty for both lepton species is related to particle identification effi- 
ciency. As the moments are ratios of measured quantities, the effects of several uncertainties, 
which are nearly independent of the of the energy, cancel out. The overall systematic un- 
certainties are 0.28% for Rg^^ and 0.06% for R^^^ for the sample, and 0.32% and 0.06% 
for the fi^ sample. These are comparable to the statistical uncertainties. Since the two 
moments are extracted from the same spectra, we must use the covariance matrix Vq,/? to 
extract the HQE parameters. Table y shows this covariance matrix for electrons and muons. 

2.4 • 10"°^ 1.1 • 10"°^ \ 

1.1 • 10-°^ 6.7-10-°^ ] 

5.2 • 10"°^ 2.2 • 10"°^ \ 
2.2 • 10'°^ 1.3 • 10'°^ ) 

1.9-10-°*^ 8.4- lO""'^ \ 
8.4 • 10"°^ 5.0 • 10-°^ j 

TABLE IL Covariance matrices for the statistical errors on Rn'^^ and Rf^^ moments. 



Electrons 



Muons 



Combined 



III. DETERMINATION OF THE HQE PARAMETERS A AND Ai 

We can determine A and Ai using the published expressions for the moments Rq and Ri 
in terms of the HQE parameters [0, ^. We use the terms describing the integral over the 
b — > qiu processes. These theoretical expressions include terms accounting for electroweak 
radiative effects and unfolding from the laboratory to the rest frame. We do not include these 
terms because our data are corrected for these effects. The theoretical expressions include 
terms through order in the non-perturbative expansion 0, ^, ^. These introduce 6 

additional form factors, denoted as pi, p2, ti, T2, t^, and t^. From dimensional arguments, 
they are expected to be of the order Aq^^ and thus they are generally assumed to be < (0.5)^ 
GeV^. In addition, pi is expected to be positive from the vacuum-saturation approximation 
p!6| . Furthermore, as Gremm and Kapustin have noted Q, the B*-B and D* — D mass 



splittings impose the constraint 



_ K{m,)Ml^MB{MD + A) - Ml^MpiniB + A) . . 

Mb + A - kK) ■ (Mfl + A) ' ^ ' 
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where n{mc) = [as{mc) /as{'m'b)]^^^^°^ and AMb^d represents the vector-pseudoscalar meson 
sphtting in the charm and beauty sectors. 

The values of the HQE parameters and their experimental uncertainties are obtained by 
calculating the from the measured moments Rq^^ and R^^^ and the covariance matrix 

a=l /3=1 



X 



) exp 



Rl^) V 



R 



th\ 



(8) 



a=0 13=0 

In Fig. ^ we show the Ax^ = 1 contours for electrons and muons corresponding to the 
quoted experimental uncertainties. 

The theoretical uncertainties on the HQE parameters are determined by varying, with 



0.09 ±0.02, 
= ± (0.5)3 



flat distributions, the input parameters within their respective errors: ly^] = 
as = 0.22 ± 0.027, A2 = 0.128 ± 0.010 GeV^, pi = ^{0.5f ± 1(0.5)3 Q^yf^ 
GeV^, Tj = 0.0 ± (0.5)3 QgY3_ 'jj-^g contour that contains 68% of the events is shown in Fig. 
^. This procedure for evaluating the theoretical uncertainty from the unknown expansion 
parameters that enter at order l/M^ is similar to that used by Gremm and Kapustin 0] 
and Bauer and Trott but different from the procedure used in our analysis of hadronic 
mass moments ||T0|- The dominant theoretical uncertainty is related to the 1/Mg terms in 
the non-perturbative expansion discussed before. 



The measured Ai and A are given in Table III 





Ai(GeV2) 


A(GeV) 


Electrons 


-0.28 ± Om\stat ± 0.06\syst ± 0-lMth 


0.41 ± OM\stat ± Om\syst ± 0.121^;, 


Muons 


-0.22 ± OM\stat ± Om\syst ± O.U\th 


0.36 ± Om\stat ± 0.08|,j,,t ± O.U\th 


Combined 


-0.25 ± om\stat ± om\syst ± o.u\th 


0.39 ± Om\stat ± OMUyst ± 0-12\th 



TABLE III: Measured Ai and A values, including statistical, systematic, and theoretical errors. 



A previous CLEO measurement used the first moment of the hadronic recoil mass [0 and 



the first photon energy moment from the 6 — 57 process . Figure ^ shows a comparison 
of our results with the previously published ones. We overlay the ellipse from the electron 
spectrum measurements extracted using \ y^\ = 0.07 for consistency with their assumptions. 
The agreement is good, although the theoretical uncertainties do not warrant a very precise 
comparison. 

Using the expression of the full semileptonic decay width given in Eq. |l], we can extract 
\Vcb\. We use F^f^ = (0.43 ± 0.01)xlO-i°MeV Q. Assuming the validity of quark-hadron 
duality, 

|V,fe| = (40.8 ± 0.5|r,, ± OA\^„a ± 0.9|t;,) x lO'^, (9) 

where the first uncertainty is from the experimental value of the semileptonic width, the 
second uncertainty is from the HQE parameters (Ai and A) and the third uncertainty is the 
theoretical uncertainty obtained as described above. 

Finally, we use the formalism of Ref. P] to extract the short range mass of the 6-quark 
mP, defined as = Mb — A^'^(GeV). Table |V| summarizes the measurement of these 



parameters for electrons, muons and the sum. The theoretical uncertainty is extracted using 
the method described above. Our result is in good agreement with a previous estimate of 
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g] derived from T(15) data, mf= 4.69 ± 0.03 GeV/c^. 
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Ai^(GeV) 


mP(GeV/c2) 


Electrons 


0.52 ± OMUtat ± 0.06|,,y,f ± O.U\th 


4.79±0.07|exp±0.11|t;, 


Muons 


0.46 ± 0.05l,tat ± 0.08\syst ± O.lllt/, 


4.85±0.09|e^p±0.11|t,, 


Combined 


0.49 ± 0.03|,tat ± 0.06|,j,,t ± O.llltfe 


4.82±0.07|exp±0.11|tfe 



TABLE IV: The measured A^*^ and ml^ . The quoted errors reflect statistical, systematic and 
theoretical uncertainties, respectively. 



IV. CONCLUSION 

Wc have measured the moments Rq and Ri of the lepton energy in ^ Xiu decays and 
obtained the HQE parameters A = 0.39±0.03|,tat ±0.06|,j,,t ±0.12|t;, GeV and Ai = -0.25 ± 
0.02|^tat ± 0.05\syst ± 0.14|t,, GeV^. These results imply that = 4.90 ± 0.08\exp ± 0.13|t/, 
GeV/c^. The short range mass ml^ is found to be 4.82 ± 0.07|exp ± 0-ll\th GeV/c^. Finally, 
we obtain |Vc6| = (40.8 ± 0.5|r,, ± 0.4|;,^ ^ ± 0.9]^^) x 10"=^. 

These results are consistent with previous measurements, and thus we have no evidence 
for quark-hadron duality violations. A quantitative assessment of these unknown errors will 
be possible only when the uncertainties are known more precisely. 
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FIG. 1: Raw electron (top) and muon (bottom) spectra from T(45) d and A ) and scaled 
continuum (o and V) data. Positive (red) and negative (green) leptons are shown separately. 
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FIG. 2: Contributions to the electron (top) and muon (bottom) spectra from background processes 
that were estimated with data. 
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FIG. 3: Contributions to the electron spectrum from background sources that were modeled with 
Monte Carlo. 
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FIG. 4: Corrected electron (green A) and muon (red D momentum spectra in the S-meson rest 
frame. 
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FIG. 5: Constraints on the HQE parameters Ai and A from our measurements of the moments 
Ro and Ri. The contours represent Ax^ = 1 for the combined statistical and systematic errors on 
the measured values. The parameters Ai and A are computed in the MS scheme to order 
and /JqcKs- 
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FIG. 6: The constraints from our combined electron and muon Rq and Ri moments, with Ax^ = 
1 contours for total experimental and theoretical uncertainties. The parameters Ai and A are 
computed in the MS scheme to order and /JqcKs- 
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compared with the combined electron and muon Rq and Ri constraints. The parameters Ai and 
A are computed in the MS scheme to order 1/M^ and Poo^- 
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FIG. 8: The combined electron and muon Rq and Ri constraints, with Ax^ = 1 contours for total 
experimental and theoretical uncertainties in the plane of Ai and A^*^, according to the scheme of 
Ref. § to order 1/M|. 
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